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ABSTRACT 



Aims. The Jk = lo-Oo rotational transition of phosphine (PH3) at 267 GHz has been tentatively identified with a Tmb ~40 mK spectral 
line observed with the IRAM 30-m telescope in the C-star envelope IRC +10216. 
Methods. A radiative transfer model has been used to fit the observed line profile. 

Results. The derived PH3 abundance relative to Ht is 6 x 10"', although it may have a large uncertainty due to the lack of knowledge 
about the spatial distribution of this species. If our identification is correct, it implies that PH3 has a similar abundance to that reported 
for HCP in this source, and that these two molecules (HCP and PH3) together take up about 5 % of phosphorus in IRC +10216. The 
abundance of PH3, as that of other hydrides in this source, is not well explained by conventional gas phase LTE and non-LTE chemical 
models, and may imply formation on grain surfaces. 
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1. Introduction 



IRC +10216 — stars: carbon — radio lines: stars — astrochemistry — line:identification — stars: 



Of the nearly 150 molecules detected so far in interstellar and 
circumstellar media, around 3/4 can be formed from just four 
elements (H, C, N, and O). The remaining 1/4 contain metals 
(Na, K, Al, Mg, and Fe), halogens (F and CI), and to a large 
extent the second-row elements Si, P, and S. The scarcity of 
molecules containing the second-row elements Si, P, and S, com- 
pared to their first-row analogues C, N, and O, reflects on one 
hand a lowe r cosmic abundance (Si/C ~ 1/8, P/N ~ 1/300, and 
S/O ~ l/30: lAsplund et alj |2005) and on the other hand impor- 
tant chemical diflferences. Second-row elements have generally 
a larger refractory character, i.e. they tend to form solid conden- 
sates and deplete f rom th e gas phase dFieldl [T974t iTurneii 1 1 99 U 
iLodders & Feglev. fT999h . Also, chemical bonds formed by first- 
row elements are generally stronger and the resulting molecules 
are more stable than the corresponding second-row analogues. 
For example, the high stability of C-C bonds allows to form a 
large variety of organic molecules, while a Si-based chemistry is 
much more limited. 

Concerning phosphorus, the very limited number of inter- 
stellar and circumstellar P-bearing molec i iles known for a long 
time, just PN and CP jTurner & BaUvl Il987t IZiurvsl Il987t 
iGuelin et all 1 19901) in contrast with the nearly 50 N-bearing 
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molecules found, may have been related to a sort of mixing of 
all the factors mentioned above. First, the more than a hundred 
fold decrease in the cosmic abundance of P compared to that of 
N. Second, the likely strong depletion of phosphorus suggested 
by observations of some P-bearing species (Turner et al., 199§. 
And lastly, some important differences between the chemistry of 
P and that of N, that are indicated by, for example, the non re- 
activity of PH;t (n = 0,3) i ons with H2 as op posed to the case of 
NH+ (iThorne et al.L IT984t lAdamsltaTl (1990). 

In recent years, however, we have seen the discoveries, 
mostly in circumstellar me dia, of new phosphorus compounds 
such as HCP PO and C CP (fAgtindez et al.', '2007"; ' Milam et all 
I2OO8I; iTenenbaum etaLl |2007;'Halfen et al., 2008), which per- 
mit to discuss the chemistry of phosphorus on a wider obser- 
vational basis. In this letter we present the tentative detection 
of PH3, the phosphorus analogue of NH3, in the carbon-rich 
circumstellar envelope IRC +10216. This species, known to be 
abundant in t he atmospheres of the giant gaseous planets Jupiter 
and Saturn ( Weisstein & Serabvnl [19961) . has never been ob- 
served outside the Solar System. 



2. Observations and results 

The phosphine molecule, PH3, is an oblate symmetric top, thus 
its rotational levels are given by two quantum numbers (J, K), 
and radiative transitions are only allowed within levels of the 
same K ladder (AJ - 1, AK - 0). The K ladders are grouped into 
two distinct forms: ortho (K - 3n, n an integer) and para (K 
3n) between which both radiative and collisional transitions are 
severely forbidden. Its rotational spectrum has been extensively 
investigated in the laboratory, allowing for the very weak "for- 
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Table 1. Observed line parameters in IRC +10216 
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SiS v=4 I 
J-14-13 'I 

1 CCC'^CH 

N = 27-26 




-30 
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Fig.l. Spectra of IRC +10216 observed with the IRAM 30-m 
telescope showing on top the Jk - Iq-Qq line of PH3 blended 
with a narrower line assigned to the J =15-14 line of SiS v = 4. 
Shown below are the immediate lower rotational transitions of 
SiS V = 4. Fits to the line profiles using the CLASS shell method 
are shown as grey thick lines. 



bidden" transitions (AJ - 0, AK - ±3) to be measured and for the 
hyperfine str ucture due to the 'H and the ^'P nuclear spins to be 
resolved (see lCazzoh & f^izzarini 2006 and references therein). 
In contrast to NH3, no evidence of inversion doubling has been 
found in the case of PH3. The electr ic dipole momen t has been 
measured as 0.57395 +0.0003 D (Da vies et al.Lll971h . 

Following the rec ent detection of HCP in IRC +10216 
dAgundez et all l2007l) . it was speculated that PH3 might be de- 
tectable in this source if the PH3/HCP abundance ratio is similar 
to the NH3/HCN one. Prompted by this hypothesis we searched 
for the fundamental ortho-PH3 Iq-Oq line, at 267 GHz, with 
the IRAM 30-m telescope. Preliminary observations with a low 
spectral resolution, 1.25 MHz, were done in 2007 May and we 
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Jk = lo-Oo 


266944.514 


266944.5(3) 


14.1(3) 


1.002(40) 


SiS V = 4 










J= 15-14 


266941.754 


266942.4(7) 


4.0(7) 


0.237(40) 


J= 14-13 


249155.372 


249155.6(4) 


3.5(4) 


0.172(40) 


J= 13-12 


231366.976 


231367.2(3) 


2.7(3) 


0.256(30) 


J= 12-11 


213576.710 


213576.9(5) 


3.5(5) 


0.207(40) 



Note.- Number in parentheses are Icr uncertainties in units of the 
last digits. Observed frequencies are given in the rest frame of IRC 
+ 10216 assuming a systemic velocity of Vlsr = -26.5 km s"'. 
" Vexp is the half width at zero level. 



found significant emission at the frequency of the PH3 Iq-Oq line 
and with the expected linewidth (Va vp - 14.5 km s~' for rn ost of 
the molecular lines in IRC +10216: lce"rnicharo et aT]|2000l) . The 
line, however, showed a profile unusual for IRC +10216: nei- 
ther really U-shaped nor flat-topped. Encouraged by this result 
we returned to the 30-m telescope in 2008 February and April to 
re-observe this transition with a 4 times higher spectral resolu- 
tion (320 kHz). Two SIS receivers operating at 1 mm were used 
simultaneously with upper side band rejections of ~ 10 dB. The 
local oscillator was shifted by 80 MHz to identify any contri- 
bution from the image sideband. An autocorrelator was used as 
backend to provide the required spectral resolution of 320 kHz. 
The secondary mirror was wobbled by +90" at a rate of 0.5 Hz. 
The pointing and focus of the telescope were checked every 1-2 
hours on Saturn, which was closer than 10° to IRC +10216. Here 
we express the intensity scale in units of main beam brightness 
temperature Tmb- The parameter Befj/Feff, used to convert T^ 
into Tmb, is 0.51 and the beam size is 9" at 267 GHz for the 
30-m. 

The resulting spectrum at 267 GHz is shown on top of Fig.[T] 
The line profile indicates the blend of a normal cusped line of 
width ~29 km s"', characteristic of optically thin lines aris- 
ing in the outer envelope (e.g. C2H), with a narrower (~6 km 
s"^) sharply peaked line, characteristic of vibrationally excited 
lines arising close to the star. The latter line was soon identi- 
fied as the J = 15-14 transition of SiS in its v - 4 vibrational 
state, bas ed on SiS laboratory spectroscopic data dSanz et all 
2003; Miiller et al.i, 2007). Emission of SiS in vibrationally ex- 
cited states up to v = 3 has been reported toward IRC +1 02 1 6 



(lTurneiill987tlFonfna et all 120061: ICernicharo et al.Ll2000h . The 
linewidths of vibrationally excited SiS are unusually narrow im- 
plying that the emission arises from the innermost envelope 
where the gas has not yet reached the terminal expansion ve- 
locity of 14.5 km s To constrain the possible contribution of 
the J = 15-14 fine of SiS v = 4 to the 267 GHz fine, we observed 
the immediate previous J transitions. In Fig. [T] we show the J = 
14-13 and J = 13-12 lines of SiS in the v = 4 state (both observed 
with a spectral resolution of 320 kHz at the time of the 267 GHz 
observations) and the J = 12-11 line (observed in a previous run 
in 2005 January with a spectral resolution of 1.25 MHz). 

Except for the J= 14-13 line of SiS v = 4, which is partially 
blended with a fine structure component of the N = 27-26 tran- 
sition of CCC'^^CH, the observations clearly show that SiS v 
- 4 lines are narrow, with expansion velocities of about 3 km 
s"' (see Table[TJ. With this in mind, the observed emission fea- 
ture at 267 GHz has been fitted by two line components, a nar- 
row one corresponding to the J = 15-14 transition of SiS v = 4 
and a wider line whose width, Vexp - 14.1 +0.3 km s"', agrees 
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with the expansion velocity of 14.5 km s ' in IRC +10216, and 
whose center rest frequency, 266944.5 +0.3 MHz, is in very 
good agreement with the laboratory frequency of the Iq-Oq line 
of PH3 (see Table [1]). The hyperfine structure is not resolved as 
the components are separate d by less than 0.2 km s ' in velocity 
jCazzoli & Puzz^in ill2006l). which is smaller than the spectral 
resolution. The good agreement between the observed and labo- 
ratory frequencies is the strongest evidence of PH3 detection in 
IRC +10216. 

Besides the SiS v = 4 line, there are some other lines with 
frequencies close to that of PH3 Iq-Oq. Most of them, such as 
SO2 309,21-3 1 8,24 at 266943.344 MHz or CH3CH2CN 154,12- 
152,13 at 266951.639 MHz, are ruled out as likely contributors 
since many other lines of these species should have been de- 
tected. A more plausible species is HC3N in the V7 = 4 vibra- 
tional sta te, whose tr ansition J = 29-28 ^ = Oe lies at 266943.313 
MHz (MboseTetaLlllOOO). Although several Hnes of HC3N in 
the vibr ational excited state vj - 1 have been observed in IRC 
+ 10216 dCernicharo et al.Ll2000h . we rule out that HC3N V7 = 4 
is the main contributor to the 267 GHz emission based on the 
upper limit of Tmb < 0.01 K that we have from our 30-m data 
archive for the lower-J transition J - 26-25 ^' = Oe transition at 
239370.171MHz. 

We tried to confirm the identification of PH3 by observing 
other transitions. The J = 2-1 line at 534 GHz is not reachable 
from the ground due to severe atmospheric absorption. We, thus, 
searched for the J - 3-2 transition at 800 GHz with the Caltech 
Submillimeter Observatory (CSO). The observations were car- 
ried out in 2008 January using the chopping secondary mode 
with a throw of ±90" at a rate of 1.2 Hz. The SIS receiver 
was tuned in double sideband and a Fast Fourier Transform 
Spectrometer was used as backend to provide a spectral reso- 
lution of 0.12 MHz. The pointing of the telescope was checked 
on Saturn. The beam size of the CSO at 800 GHz is 11.5" and 
the beam efficiency is 0.28. In spite of the good atmospheric con- 
ditions (zenith sky opacity at 225 GHz was 0.04-0.08), the high 
opacity of the atmosphere at 800 GHz (Tsys ranged from 3000 K 
to 8000 K) did not allowed to reach a noise level low enough to 
confirm or discard the presence of PH3 (see top panel in Fig.|2]i. 

3. Analysis and discussion 

To interpret the observations we have computed line pro- 
files by means of excitation and radiative transfer calcula- 
tions based in the Large Velocity Gradient (LVG) formalism 
( CastoA 19701), co upled to the spectral catalog of J. Cernicharo 
dCernicharo et al.L 2000). We consider separately both the ortho 
(0-PH3) and para (P-PH3) species of phosphine (the ortho-to- 
para ratio is assumed to be 1, the statistical value). The energy 
levels and transition frequ encies are computed from t he rota- 
tional constants reported bv lCazzoli & Puzzarinil d200 6l). We in- 
clude rotational levels in the ground vibrational state up to Jk 
= 76 for 0-PH3 and up to ]k = 55 for P-PH3. As rate coeffi- 
cients for collisional de-excitation of o-PH3(p-PH3) with H2 and 
He, we have ado pted those compute d for col Usions of o-NH3(p- 
NH3) with P-H2 dPanbvet al.lll988h and He ( Mach in & Rouefll 
l2005h respectively, properly corrected to the case in which in- 
version doubling is not resolved. 

We assume a distance to IRC +10216 of 150 pc and sim- 
ulate the circumstellar envelope as a spherically distributed 
gas expanding at a constant velocity of 14.5 km s"^ The 
gas density and tempe rature radial profiles are taken from 
lAgundez & Cernicharol(l2006 ). PH3 is assumed to be distributed 
between an inner radius Rin and an outer radius Rout, with a con- 
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Fig. 2. The lower panel shows the PH3 Jk = Iq-Oq line at 267 
GHz as observed with the IRAM 30-m telescope, in which the 
fit to the SiS V = 4 J = 15-14 line has been subtracted. The upper 
panel shows the observation at 800 GHz with the CSO smoothed 
to a spectral resolution of 2 MHz. The grey thick lines in both 
panels correspond to the line profiles given by the LVG model. 



stant abundance, Jt:(PH3), relative to H2. We take Rjn = 1.3 x 10^^ 
cm (about 20 stellar radii), the value adopted by lHasegawa et al] 
(2006) for NH3. The values of Rout and x(PH3) are varied un- 
til obtaining the best fit to the observed PH3 Iq-Oq line profile. 
Since both gas density and temperature vary greatly with radius, 
we have divided the envelope into various shells and solved for 
the level populations in each shell independently of the others. 
We then compute the emergent intensity and weight it with the 
main beam of the selected telescope. 

In the bottom panel of Fig. |2] we plot the observed Iq-Oq 
PH3 line profile, obtained by subtracting the fit to the SiS v 
- 4 blended line. Also shown is the profile resulting from our 
best LVG model, whose parameters are Rout = 2 x lO"' cm and 
ji[:(PH3) = 6 X 10 The model indicates that the upper level in- 
volved in the Iq-Oq PH3 line is mostly populated by collisions in 
the inner 4 x lO'^ cm, which corresponds to an angular diameter 
of about 4" and is thus spatially diluted in the 9" beam of the 
30-m telescope. The top panel of Fig. |2] shows the correspond- 
ing prediction for the PH3 J = 3-2 line profiles as observed with 
the CSO. It is seen that the expected intensity is well below the 
noise of the observed spectrum and thus is consistent with the 
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IRAM 30-m result for the J = 1-0 line. The spatial dilution is 
also important for the CSO observations and plays against our 
efforts to detect the J = 3-2 lines. 

It should be noted that the abundance and distribution of PH3 
derived are very uncertain. Since we have just one single PH3 
line observed, there exists a considerable degeneracy between 
models with different values of Rin, Rout and ^(PHs). Moreover, 
infrared pumping to excited vibrational states, not considered in 
our model, may play an important role in the excitation of the 
rotational levels in the ground vibrational state. PH3 has indeed 
many vibrational bands in the spectral region around 10 fim, a 
wavelength at which th e central source has its maximum flux 
dCemicharo et al.Lll999l) . 

Confirmation of our tentative detection may rely on further 
observations of other PH3 transitions. The J = 2-1 and J= 3- 
2 lines, at 534 GHz and 800 GHz respectively are observable 
with the Herschel Space Observatory (HSO), although our LVG 
model predicts intensities somewhat weaker than with the CSO 
due to the larger dilution effect. In the case of the Atacama Large 
Millimeter Array (ALMA) the high angular resolution that it will 
provide, better than 3" at 800 GHz, will perfectly fit with the 
expected size of the J = 3-2 emission. The predictions indicate 
main beam brightness temperatures of about 1 K, which pro- 
vide a good opportunity of detection in spite of the high atmo- 
spheric opacity at this frequency. We note that these predictions 
are based on our best LVG model, but the J = 2-1 and J - 3-2 
lines could be more intense than expected if infrared pumping is 
playing an important role, or if PH3 is present in hotter regions, 
inner than 1.3 x 10'^ cm. 

The derived abundance of PH3, x(PH3) = 6 x 10"^, implies 
that it takes up about 1.3 % of the avai lable phosphorus. Th e 
PH3/HCP ratio is then found to be 1/2.3 dAgundez etaUl2007[) . 
noticeably hig her than the NH 3 /HCN ratio, that is about 1 /50 in 
IRC H-10216 (iHasegawa et all 120061 iFonfna et all l2007h . The 
formation of PH3 in the gas phase is difficult to explain. Its 
LTE abundance in the inner envelope is rather low, < 10 and 
gas p hase reactions yield n o net formation in the outer enve- 
lope dAgundez et al.L |2007[) . Besides PH3, other hydrides such 
as NH3, CH4, H2O, SiH4, and H2S a re observed in IRC -H i 02 16 
with relatively high abundances ([K eadv & Ridg wavi 119931; 
IHasegawa et all 120061: lAgundez & Cernicharo, 20061) . Many of 
them are usually assumed as parent molecules, i.e. formed in 
the inner en velope, in most chemical m odels of IR C +10216 
(e.g. MacKa v & Charnle v 2001; Agundez & Cernic haro 2006). 
However, c hemical equilibrium c alculations, similar to those 
reported in lAgundez et al.l d2007l) . indicate abundances much 
lower than observed, except perhaps for CH4. A widely invoked 
explanation, when gas phase chemistry fails to explain an ob- 
served abundance, is that of grain surface reactions. In the case 
of hydrides such as PH3, a likely formation process is the direct 
hydrogenation of the heavy atom taking place on grain surfaces. 

4. Conclusions 

We have tentatively detected PH3 in IRC +10216 through its 
Jk = Iq-Oo transition at 267 GHz. The derived abundance rela- 
tive to H2 is 6 X 10"**. Although with a considerable uncertainty, 
this value is similar to the HCP abundance found in this source 
dAgundez etaUl2007l) . These two species, HCP and PH3, would 
then take up about 5 % of the phosphorus in IRC +10216. The 
formation of PH3, unlike that of HCP, is difficult to explain in the 
gas phase and could occur on grain surfaces. It remains a target 
for the future to confirm this tentative detection by observing the 
J = 3-2 transition, at 800 GHz, with the ALMA facility. Also, 



further observations at 267 GHz in other s ources such as CRL 
2688, where HCP has been also detected dMilam et all |2008|) . 
will be of great interest to support this tentative detection and to 
understand the chemistry behind it. 
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